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5 CROSS-REFERENCE TO RELATED APPLICATION 
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7 This application claims the benefit of U.S. provisional patent application Serial No. 

8 60/270,960 filed February 22, 2001 . 
9 

10 BACKGROUND OF THE INVENTION 
11 

12 The present invention relates to magnetic resonance imaging and, in particular, to 

1 3 gradient coils used therewith. 

, 1 4 A Magnetic Resonance Scanner where the direction of the main static field is orthogonal 

Q5 to the pole surface and to the magnet's main symmetry plane is commonly called an Open 

;;§6 Magnet System. 

.: 

: : 

£17 Previous methods for production of linear magnetic gradients in such MRI systems 



Ug8 consist of winding discrete coils, in a bunched or distributed fashion on an electrically insulating 

„ 19 hollow right cylinder coil-form, and driving said coil with a limited voltage. Conventional 

!;;|0 bunched coil designs are the Maxwell and Modified Maxwell Pair for the Z-gradient coil 

j'igl (parallel to the main magnetic filed), and the Golay or Modified Golay (multi-arc) Saddle coils 

!=-J2 for the X/Y gradient coils. Previous methods consisted of iteratively placing elemental loops or 

1 ! £3 arcs on the cylindrical coil form until the desired gradient strength, gradient uniformity and 

24 inductance were achieved. In general, these types of designs constituted as method referred as a 

25 "forward approach" whereby a set of initial coil positions (i.e. initial current distribution) were 

26 defined. The gradient field, the energy/inductance were calculated and if they did not fall within 

27 specified design criteria, the positions of the coils will be shifted (statistically or otherwise) and 

28 the new outcome will be re-evaluated. This iterative procedure will continue until a suitable 

29 solution that satisfies all the design criteria is found. 

30 An alternative approach in the design of gradient coils is referred as an "inverse 

3 1 approach". According to this methodology, a set of predetermined constraint points for the 

32 gradient magnetic field is chosen inside the desired imaging volume. Via a Lagrange 

33 minimization technique involving the stored magnetic energy or power of the structure, the 

1 



1 appropriate continuous current distribution that generates the predetermined gradient field inside 

2 the imaging volume is derived. Using a stream function discretization technique, the continuous 

3 current distribution of the coil can be approximated by a set of discrete current loops that are 

4 connected in series and share the same constant current. In order to ensure that the discretization 

5 process is valid, the gradient field is regenerated using the Biot-Savart law to the discrete current 

6 distribution and compare the results with the ones of the analytic approach. If the error is 

7 relatively small, the discrete current distribution is accepted as a valid one. 

8 Open magnet geometries with vertically directed fields have been a successful presence 

9 in the MRI arena due to the inherent open architecture scheme that promotes patient comfort and 

10 diminishes patient claustrophobia. Recently, vertical field systems with higher main field 

1 1 strengths (0.5 T to 1.0 T) were also introduced in an attempt to improve image quality, as well as 
: 12 to expand the spectrum of clinical applications, especially fast imaging techniques not available 
i;D3 on lower field strength MR scanners. The necessity for utilizing such imaging techniques is 

: : |4 strongly dependent upon the availability of a suitable gradient system capable of delivering high 

;; : ?5 gradient field strengths with a significant improvement to the gradient coil's slew rate. For open 

\ M6 magnet geometries with vertically directed fields, biplanar gradient coils have been the geometry 

7 of choice, since there are conformed to the main magnet's geometric shape and assist in 

;;l8 improving patient comfort and reducing patient claustrophobia. There is however a drawback in 

n39 regards to the biplanar designs. Because of the demand that the coil must have enough gap to 

: i|0 allow for clear access to the human torso, the distance between the two planes is excessively 

fill high. Increasing the gap between the two planes reduces the efficiency of the biplanar gradient in 

22 a level. Thus for a whole body gradient coil set, it is very difficult to generate gradient fields of 

23 the order of 50-60 mT/m with slew rates approaching 1000 T/m/sec. An alternative solution to 

24 this issue is the introduction of uniplanar gradient coil designs, which offer superb gradient field 

25 characteristics over a very limited imaging volume. Uniplanar gradients are mainly utilized on 

26 local imaging applications and are not recommended for applications with sizable FOV such as 

27 the one covering the human head. One viable alternative for performing ultra fast imaging of the 

28 human head for a vertical field system is the use of insertable biplanar gradients in a very close 

29 proximity to the human head. Due to the reduction of the gap between the planes the efficiency 

30 of the coil is significantly improved, and such a design is capable of generating high peak and 

31 high uniformity gradient fields over a volume similar to the size of the human head (25 cm 
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1 DSV). But such a design is prohibitive because it violates the two principles of which the open 

2 magnet designs are based on. Placing the two planes over the human head and especially over 

3 the patient's eyes, it compromises the patient's comfort and enhances patient anxiety. 

4 Gradient coil sets in common use make serious tradeoffs between coil performance and 

5 maintaining the "open," non-claustrophobic configuration, desired in open MRI systems. 
6 

7 SUMMARY OF THE INVENTION 

8 

9 A MRI gradient coil set includes a uniplanar Z-gradient coil, a biplanar X-gradient coil, 

10 and a biplanar Y- gradient coil. The coil set provides an open Z-axis face. 
11 

12 BRIEF DESCRIPTION OF THE DRAWINGS 

;*? 

m 



;A FIG. 1 is a perspective view of a possible gradient coil set according to the invention. 

i-EB FIG. 2a is a rable showing on-axis linearity, 

i jjS FIG. 2b is a table showing off-axis uniformity. 

?7 FIG. 3 is a table showing on-axis field behavior. 

i;U8 FIG. 4 is a table showing on-axis field behavior, 

j'jft? FIG. 5 is a gradient pattern. 

j §0 FIG. 6 is a gradient pattern, 

i |1 FIG. 7a is a gradient pattern. 

22 FIG. 7b is a gradient pattern. 

23 FIG. 8 is a table snowing on-axis field behavior. 

24 FIG. 9 is a gradient pattern. 

25 FIG. 10 is a gradient pattern. 

26 FIG. 11 is a gradient pattern. 

27 FIG. 12 is a gradientlpattern. 

28 FIG. 13 is a gradient pattern. 

29 FIG. 14 is a table showing on-axis behavior. 

30 FIG. 15 is a gradient pattern. 

3 1 FIG. 1 6a is a gradient pattern. 



3 




1 FIG. )i 6b is a gradient pattern. 

2 FIG. l\a is a gradient pattern. 
3 . FIG. 1 7ft is a gradient pattern. 

4 FIG. 18 i^a perspective view of a possible gradient coil set according to the invention 

5 with shielded biplamar coils. 
FIG. 19 is a perspective view of a possible gradient coil set according to the invention 

with unshielded biplaimr coils. 

FIG. 20 is a perspective view of a possible gradient coil set according to the invention 

9 with shielded biplanar coils and incorporating radio frequency coils. 

10 FIG. 21 is a perspective view of a possible gradient coil set according to the invention 

1 1 with unshielded biplanar coil&and incorporating radio frequency coils. 

12 FIG. 22 is a perspectiveYiew of a possible gradient coil set according to the invention in 

:: t$> 

i;g3 an array configuration. 

i'ljl FIG. 23 is a perspective vieV of a possible gradient coil set according to the invention 

; ! ;» incorporating knee and upper thigh fl^equency coils. 

; if|6 FIG. 24 is a perspective view of a possible gradient coil set according to the invention 

n \ 

T7 incorporating knee and foot radio frequency coils. 

Hi 

! |9 DESCRIPTION OF THE PREFFERRED EMBODIMENTS 

So 

!::f! 

f il In this invention, a new alternative design for an insertable open architecture gradient coil 

22 set is given which is focused on enhancing patient comfort and abolishing patient anxiety, while 

23 delivering a gradient field with high peak strength and improved slew rate. According to this 

24 architecture, the Z gradient coil is a uniplanar design that is placed in parallel with the main 

25 magnet's poles. In order to form a complete gradient set, the two transverse (X, Y) gradients are 

26 generated by a pair of biplanar designs which are placed perpendicular to the magnets poles, but 

27 they may be physically mounted on the same side with the uniplanar Z gradient. The gradient 

28 field generated from the two transverse coils exhibits similar peak and slew rate performances as 

29 the uniplanar Z gradient coil over a comparable imaging volume. When this coil is designated for 

30 head imaging applications, parabolic cut-offs on the two side transverse biplanar coils may be 

31 implemented in a symmetric fashion. Such an insertable coil design does not disrupt the vision of 
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1 the human subject and thus makes it more patient friendly without imposing restriction on the 

2 field strength and quality. The proposed design can also be implemented in a self-shielded 

3 configuration for reduction of induced eddy current effects on the vicinity of the main magnet's 

4 poles. In addition, the proposed design can be paired with an additional set of gradient coils 

5 (shielded or not) in an array configuration with a capability of extending the imaging coverage of 

6 the gradient coil system, while improving patient comfort and diminishing patient 

7 claustrophobia. 
8 

9 Theory 

10 

1 1 FIG. 1, shows the geometrical shape for a self-shielded open architecture gradient set. 

12 The gap between the two primary planes on the biplanar section of the coil is denoted as 2a while 

l. .JL 

l§ the distance between the shielding planes of the biplanar section of the coil is denoted as 2b. For 

i:# the uniplanar section of the coil, the distance between the plane section of the coil and the 

i'l geometric center of the coil is denoted as a', while the distance between the geometric center of 

"4 

M the coil and shielding part of the biplanar section as b\ The schematic representation of a head 

rl 

17 Transmit/Receive Radiofrequency coil is also included. 

StP Axial (Z) Uniplanar Gradient 

I 

M Let us assume that the position of the plane for the primary gradient coil is at z=-a' from 

22 the geometric center of the coil, while the position of the shielding coil is z=-b'from the center of 

23 the coil. The current density distribution for this design is confined on the xy plane and has the 

24 expression: 

Ax,y) = [J x (x,y)x + J Y (x,y)y] (1) 

25 
26 
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1 In addition the expression for the vector potential is given by: 
2 



4 -A. f r^_L^^L^(^^^^^ (2) 

4 =0" j J>t//?-pL^ (3) 

3 

4 And for the region outside the planes z<-b', the expression for the vector potential is: 

4 = A. J (4) 

A f Hfaff ; 1 ^r^^V/^^+e^^V/te^l/orz^^ (5) 

■jl where J( X>y ) (a,P) represent the double Fourier Transform of J( X , y ) (x,y), respectively. Since there 
! ; : |> is no charge free in the vicinity of the gradient coils, the application of the current continuity 
>',j equation (V*J=0) leads to a relationship between the components of the current density as: 

B J y {a,P) = ~J x (a,P) (6) 



il 



!;^9 From equations (2)-(5) and equation (6), the expression for the vector potential at the location of 
i iiO the two planes is given by: 

■ <af 

VU Vdadp - 1 ^ Mj;\a,P)^e-^ b '- a y:\a,P)}forz = -a (7) 

&r J ^ Ja 2 +/3 Z L J 

4 =^J J>#-jP»j=^ (8) 



4 =0 { C ^-^^^[^^(^^^/(o;/?)]^-^ 
4 = ^ I £^-^^^^ 



(9) 
(10) 



11 

12 The expression for the magnetic field is: 
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Bz =- d 4 f r dad p V*' + >*' ^ m^^^j-a- (g> p) 

^-e- 2 ^^y o rz>-a' (11) 

1 

2 

3 and since the magnetic field is symmetric along the x and y directions and monotonically 

4 increased along the z direction, equation (11) becomes: 



B z =-^\ ft/or#^^co^ (12) 
with 

j;\a,/J) = ^^ 2(b '- a) j; a \a,p)(shielding condition B z =0 for z>-b') (13) 
q and 

j;| y/(a,^)=J/(a, > 3)=>J/(a,^)=4j j°d^cos^)smO^/*W) (14) 
UU6 Furthermore, the expression of the stored magnetic energy is: 

| ^•■ 7= T 6 ^ J 7^^^^U/(«^)| 2 [l-e-^^] (15) 

nU7 In order to find the optimum continuous current density, which will generate a gradient magnetic 

ny 

rjS field in predetermined locations inside the imaging volume, can be extracted via a Lagrange 

: y 9 optimization technique, According to this technique, a functional expression, which is a 

10 quadratic function of the current density is constructed. An example of such an expression 

1 1 includes which includes the stored magnetic energy and the magnetic field, and can have the 

12 form: 

F{j;°') = W a -2>> Z (0)-2?„ C (?,0] (16) 

7=1 

13 where X y are the Lagrange multipliers and B zsc (rj) are the constrained values of the magnetic field 

14 inside the desired imaging volume. Minimizing F with respect to J x ~ a , a linear matrix equation 

1 5 for J x " a is obtained : 
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= -0Zlj cos(ca y )cos(^)^^ (2j+0,) 07) 
1 and the Lagrange multipliers are determined via the constrained equation : 



= B zSCi => C, =fV J Jjo^cosCaxJcosO^^-^^ 



(18) 



2 

3 For the case of a non-shielded uniplanar axial gradient set, the above expressions are the 

4 same at the limit when b'-»oo. 

5 

U 

i"J> From the continuous current distribution and using employing the stream function technique, a 

discrete set of loops that closely approximates the continuous current density pattern is 
£3 generated. In order to verify the accuracy of the discretization, the gradient field of the coil is re- 

'is! 

; J9 evaluated by applying the Biot-Savart law to the discrete current distribution. 

i: Po 
ill 

r §2 Transverse X or Y Uniplanar Gradient 



m 

a 

f^4 The expression for the magnetic field is from equation (1 1): 



B z =-d z A y =A J fcdad ^" 2 c im ^e^^j/{a,P) 
\\-e 2 ^^forz>-a' 

15 
16 

17 and since the magnetic field is symmetric along the x and z directions and antisymmetric along 

18 the y direction, equation ( 1 1 ) becomes : 
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^\ [jcdp^-^s\^ (12a) 



J/ 6 ' (a, P) = -e^^^j-"' (a, ^(shielding condition : B z = Ofor z > -b') (1 3a) 
and 

J x ~ a \a,p) = 4 J |°ax dy sin(ccc) sin($>) -V' (*,;>) (14a) 

1 

2 Furthermore, the expression of the stored magnetic energy is: 



W„-±Jd>x A.J=^\Jdad/3^*^\j/ (15a) 

3 In order to find the optimum continuous current density, which will generate a gradient magnetic 
f$ field in predetermined locations inside the imaging volume, can be extracted via a Lagrange 

optimization technique, According to this technique, a functional expression, which is a 
; ^6 quadratic function of the current density is constructed. An example of such an expression 
Liy7 includes which includes the stored magnetic energy and the magnetic field, and can have the 
Q 8 form: 



F{j; a ) = W m -Z^k(r y .)-5 25C (F.)] (16) 



j .«j9 where A,j are the Lagrange multipliers and B zsc (rj) are the constrained values of the magnetic field 
1 '^0 inside the desired imaging volume. Minimizing F with respect to J x " a , a linear matrix equation 

1 1 for J x " a is obtained: 

= -P&j sin(ax y )cos(^,) e -^ ( ^ a,) (17a) 

12 and the Lagrange multipliers are determined via the constrained equation : 



1% = B zsq =>C,=^VJ Idadfismiaxjcosifrje-^^ 
g^.si^orpcos^^-^'^^^l-e- 2 ^^^] (18a) 
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1 For the case of a non-shielded uniplanar axial gradient set, the above expressions are the 

2 same at the limit when b'->oo. 

3 

4 From the continuous current distribution and using employing the stream function technique, a 

5 discrete set of loops that closely approximates the continuous current density pattern is 

6 generated. In order to verify the accuracy of the discretization, the gradient field of the coil is re- 

7 evaluated by applying the Biot-Savart law to the discrete current distribution. 
8 

9 Phased Array considerations: Mutual energy between two coils. 
10 

1 1 Let us assume a second planar set with planes locations of z=-a ' and z=-b 9 . If the second 

set is shifted along the x or y directions or both by an amount of xo and yo> respectively, the 

.;£3 analytical expression of the current density distribution is given by: 

m 

'is jf'- 1 ' (a, fi) = e ip » e iax * f f dx dy e***-'* e*^** J x '' a ' A ' (x, y) (1 9) 

My 00 

ii 

t'37 Considering for simplicity of the evaluations, that the second set is shifted only along the y 

nj8 direction by an amount of yo. In this situation the mutual energy between these two modules is 

l 39 given by: 

mo 

= I \d^ x . =JL \ *]da dp J a2 +fc -^jr\a,P)Jr\a,P) 

21 V -co r 



-yla 2 + fi 2 (a'-a>) ^ _ ^a 2 */} 2 {b'-{a'-a<)) ^ ^Q) 



22 

23 If both planes represent the same gradient coil axis, equation (20) becomes 



+qo / 2 d 2 

= -^f \da dp V " + P cos(Py 0 )j; a ' {a, p)j; s " (a, p) 
24 471 i P 



-Ja 2 + /3 2 («*-a') J*j _ e -2ja 2 +/3 2 (*•-(«'-«•» J 



(21) 
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The optimum phased array configuration occurs at the yo points where the value of the mutual 
inductance between the two sets is minimized. 



Transverse (X, Y) Biplanar Gradient Coils 



Let us assume that for the self-shielded biplanar design, the gap between the two planes of the 
primary gradient coil is y=2a, while the gap between the planes of the shielding coil is z=2b. 
The current density distribution for either the primary gradient coil or the secondary gradient coil 

j(x 9 y ) = [/;• (jt, z)x + r 2 a (jc, z)z] S(y - a) + [ J~ x a (x, z)x + J~ a (x, z)z] 8{y + d) + 

[y; 6 (jc, z)Jc + y 2 + * z)f] ^(^ - b) + [/;* z) * + z)z] <?( ^ + 6) (22) 

is confined on the xz plane. The analytical expression for the current density of the coil is: 
In addition the expression for the vector potential is given by: 

where J (XfZ ) (a,P) represent the double Fourier Transform of J (X)Z )(x,z), respectively. Since there 



Mo 



&r 2 J ^ V^+A 2 



&T 2 J *• L? A 



for\y\<a (23) 



J°{a,P)+e 



Jj7ff(b+y) , -b 



(ccfi) 



for\y\>b (24) 



is no charge free in the vicinity of the gradient coils, the application of the current continuity 
equation (V»J=0) leads to a relationship between the components of the current density as: 



a 



7 2 (a,/?) = --/>,/?) 



(25) 



From equations (23)-(24) and equation (25), the expression for the vector potential at the 
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location of the two planes is given by: 



A 



fory=a 



a 



4 : 
A- 



■Jk. 



f Vdod p : 1 



^ ( ^V/(a^+ e ^ ( ^(a,/5) 



:4 f IW^LW-* 



(26) 



/orv = a (27) 



/orv=-a (28) 



1 



=4 r r dC dp r 



fory=-a (29) 



/ory=Z> (30) 



J s b <p l .P>+J*+™J x +{a t fi> 



fory=b (31) 



i4. = 



^0 

8* 2 



V« 2 +/? 2 



Jw^j; (a, + e -lw«+°)j-° (a> ^ + 
^< 2i >y>,/?)+.//V/*) 



8^ 2jJ - V« J + /? 2 l/J[e W(24, y/(a^) + y.- fr («^) 



/orj = -Z, (32) 

fory = -b (33) 



The expression for the magnetic field is 

B z =-8 2 A y — ^ J ^Jadpe iea ^ 



for\y\<a (34) 



>"M>Z> (35) 



Transverse Y Gradient Coil 
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For the Y gradient coil, the gradient magnetic field is symmetric along the x and z directions 
while it is asymmetric along the y direction. In this case the relationship between the current 
densities on the adjacent planes is: 
J x °{a,B) = j;"(a,P) and J x (a, 0) = J ~ b (a,fi) 4 

The expression of the magnetic field from equations (34) and (35) becomes: 



B=- 



%n 2 



\ £dad/3e iac+iJi 



B = 



~ d ^= S, 2 



for\y\<a (36) 



for\y\>b (37) 



The gradient field shielding condition in the region outside both planes demands that B z =0 for 
|y|>b. Equation (37) leads to: 



J x - b \a,P) = -e-^ {b '- tt) j; a \a,p) (38) 
whichleadsto: 

B z =--0- 1 £dadpcos(a>c) cosO&)e-^° sinhtyja 2 + J3 2 )J x "{cc,p^\ _ e -W^V<o j (39) 
and 

J x {a,P) = 4j %dxdz cos(ac) cos(0y) J x "(x,z) (40) 
Employing the current density's symmetry conditions, the shielding condition and the continuity 



V — oo ^ 

equation, the expression for the stored magnetic energy is given by: 



(41) 



In order to find the optimum continuous current density, which will generate a gradient magnetic 
field in predetermined locations inside the imaging volume, can be extracted via a Lagrange 
optimization technique, According to this technique, a functional expression, which is a 
quadratic function of the current density is constructed. An example of such an expression 



13 



1 includes which includes the stored magnetic energy and the magnetic field, and can have the 

2 form: 

FV x a ) = W m -XA 7 k(r y .)-^ c (0)] (42) 

3 where Xj are the Lagrange multipliers and B zsc (rj) are the constrained values of the magnetic field 

4 inside the desired imaging volume. Minimizing F with respect to J x a , a linear matrix equation for 

5 J x a is obtained: 



(.alt 



j; = - , f 2 6 . IX cos(ot ; )cos(A.)sinh(^Va 2 +p 2 ) (43) 

6 and the Lagrange multipliers are determined via the constrained equation : 



\j C,=^-J J</arf/?cos( ox,)cos( /?z,)sinh( j,,^ + P 2 ) ; f ■ , = (44) 

; y * V a + fi cosh ( * yfa^ + J*) 

rl 



sinh(( ft - a ) -Ja 2 + P 2 ) £ /I , cos( a jc y ) cos( fiz , ) sinh( y } J a 2 + P 2 ) 



if 

: ^ And the expression for the stored energy is: 



w 



8 2 -b-Ja' + p-* 

= £JL f \dadp— i =M== . sinh(( b - a)Ja 2 + fi 2 ) ( 45 \ 

m 2 J •» V« 2 + p 2 cosh( ayja 2 + P 2 ) V ( i) 

^Ji, A j cos( ajc ; ) cos( Pz } ) sinh( y V« 2 + /? 2 ) j 

11 

12 

13 For the case of a non-shielded biplanar transverse gradient set, the above expressions are 

14 the same at the limit when b->oo. 
15 
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From the continuous current distribution and using employing the stream function technique, a 
discrete set of loops that closely approximates the continuous current density pattern is 
generated. In order to verify the accuracy of the discretization, the gradient field of the coil is re- 
evaluated by applying the Biot-Savart law to the discrete current distribution. 



Transverse X Gradient Coil 



For the X gradient coil, the gradient magnetic field is symmetric along the y and z directions 
while it is asymmetric along the x direction. In this case the relationship between the current 
densities on the adjacent planes is: 

J x a {a,/5) = -j; Q {a 9 P) and J x \a,P) = -j; b {a,P) 

The expression of the magnetic field from equations (34) and (35) becomes: 



for\y\<a (46) 



B. =■ 



<ax+ipz 



J x a {a,p){-e^ a 2cos\i^a 1 + ))- 
J x \a,/3)(-e-"^ 2 2cosh^« 2 +^ 2 )) 



for\y\>b (47) 



The gradient field shielding condition in the region outside both planes demands that B 2 =0 for 
|y|>b. Equation (47) leads to: 



j;\a,P) = -e^ (b - a,) j; a \a,P) (48) 
whichleadsto : 

B z J £dadj3sm(ax) cos(/feK^ fl coshfr^ 2 + P'V^a^-e 2 ^^ (49) 

and 

J x a (a 9 /3) = LJ x a (a 9 J3) =>J x a (a,P) = 4$ [dxdz sm(ax) cos(/fc) J x \x,z) (50) 
Employing the current density's symmetry conditions, the shielding condition and the continuity 



W u =\\d>x ^*J=^ / 7^^^^l^ a («»^)! 2 [l-e-^)[l-e- 2 ^^>] (51) 
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1 equation, the expression for the stored magnetic energy is given by: 
2 

3 In order to find the optimum continuous current density, which will generate a gradient magnetic 

4 field in predetermined locations inside the imaging volume, can be extracted via a Lagrange 

5 optimization technique, According to this technique, a functional expression, which is a 

6 quadratic function of the current density is constructed. An example of such an expression 

7 includes which includes the stored magnetic energy and the magnetic field, and can have the 

8 form: 



Id 



rt3 

m 
m 
m 
a 

■ ! t? 



FV°) = W m -^AjlB^-B^j)] (42) 



7 = » 



9 where are the Lagrange multipliers and B zsc (rj) are the constrained values of the magnetic field 

.10 inside the desired imaging volume. Minimizing F with respect to J x a , a linear matrix equation for 

O J x a is obtained: 

.; k=l 

: : 

m . 

O B 2 e -"J<* 1+ 0 2 n 

J x a = r 4— ~-j—-^Aj sm(aKj)cos(/kj)cosh(yjJa 2 + 0 2 ) (52) 

fi2 and the Lagrange multipliers are determined via the constrained equation : 



r _ Mo f Idccdfi sin( ax t ) cos( /fe f ) cosh( y i : ^]a 2 + ft 2 ) —t=£= 

ij K 2 J J - • ■ ^Ct 2 +P 2 ^ 

— — ■ ^A y sin(ax 7 .)cos(^ 7 )cosh(3; y .Va 2 + fi 2 ) 

(1-e aV0f +p ) 7=i 

14 

1 5 For the case of a non-shielded biplanar transverse gradient set, the above expressions are 

16 the same at the limit when b->oo. 
17 

1 8 From the continuous current distribution and using employing the stream function technique, a 

1 9 discrete set of loops that closely approximates the continuous current density pattern is 



2 
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generated. In order to verify the accuracy of the discretization, the gradient field of the coil is re- 
evaluated by applying the Biot-Savart law to the discrete current distribution. 



Transverse Z Gradient Coil 



For the Z gradient coil, the gradient magnetic field is symmetric along the y and x directions 
while it is asymmetric along the z direction. In this case the relationship between the current 
densities on the adjacent planes is: 

J x a {a,P) = -j;\a,P) and J x \a 9 fi) = -j;\a,p) (54) 

The expression of the magnetic field from equations (34) and (35) becomes: 



B =- 



Mo 
%n 2 



J £dadpe l 



icx+ipk 



\ fcdcdpe***"* 



J x a (a,P)(e-^ (a - y) +e^ {a+y) )+ 

J x ' (a,P){-e^°2cos^cc 2 +p 2 ))- 
j/(a,/?)(V^2cosh^Va 2 +P 2 )) 



for\y\<a (55) 



for\y\>b (56) 



The gradient field shielding condition in the region outside both planes demands that B z =0 for 
|y|>b. Equation (56) leads to: 



j; b \a,P)^-e^ b, - 0,) j; a \a,fi) 
whichleadsto: 



(57) 



B * =J ^ \ £dcd/3cos(ca:) sin(fr)e-^ a cosh^a 2 +p 2 )J x a (a,p^\-e- 2 ^ (b - a) ^ (58) 



and 



J x °(a,P) = iJ x a {a,P)^>J x a {a,P) = -4i\ |° dx dz cos(ck) s\n{pz) J x a (x,z) (59) 
Employing the current density's symmetry conditions, the shielding condition and the continuity 



equation, the expression for the stored magnetic energy is given by: 
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1 In order to find the optimum continuous current density, which will generate a gradient magnetic 

2 field in predetermined locations inside the imaging volume, can be extracted via a Lagrange 

nJ/) = W m -f d A J [BA?j)-B 25C (r J )] (42) 

3 optimization technique, According to this technique, a functional expression, which is a 

4 quadratic function of the current density is constructed. An example of such an expression 

5 includes which includes the stored magnetic energy and the magnetic field, and can have the 

6 form: 

7 where are the Lagrange multipliers and B zsc (rj) are the constrained values of the magnetic field 

8 inside the desired imaging volume. Minimizing F with respect to J x a , a linear matrix equation for 

9 J x a is obtained: 



i'iji 



, f , 6 f?!,, IX- cos(ax.)sin(A 7 )cosh(^Va 2 +/? 2 ) (60.) 



-fl 

j; = " 



; i§ and the Lagrange multipliers are determined via the constrained equation : 



C _4J J^y?cos(ax.)sin(>gzJcosh(^Va 2 +yg 2 ) ^ = = ( 61 ) 



3 L L " ^A.cos(ox.)sin( > &;)cosh(>;.V« 2 +/? 2 ) 



11 

12 

13 For the case of a non-shielded biplanar transverse gradient set, the above expressions are 

14 the same at the limit when b->oo. 

15 

16 From the continuous current distribution and using employing the stream function technique, a 

17 discrete set of loops that closely approximates the continuous current density pattern is 

18 generated. In order to verify the accuracy of the discretization, the gradient field of the coil is re- 

19 evaluated by applying the Biot-Savart law to the discrete current distribution. 
20 
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1 Mutual Inductance evaluations for Biplanar Open Architecture coils 

2 

3 Let us assume a second biplanar gradient set with plane gaps for the primary and secondary coils 

4 at locations of z=2a and z=2 b , respectively . If the second set is shifted along the x or z 

5 directions or both by an amount of xo and zo, respectively, the analytical expression of the current 

6 density distribution is given by: 
7 

8 J/ 5 ' - * (a,/?) = e-'* 0 e- to ° J £<£t dze^^e"^ j;' a '~ l {x,z) (62) 
9 

10 Considering for simplicity of the evaluations, that the second set is shifted only along the z 

1 1 direction by an amount of z 0 . In this situation the mutual energy between these two modules is 

12 given by: 



13 



no 



W u ^ =1 \d'x A™ • J (5 '*> =-^- f Ida dp ^\^ 
2 1 An 1 J _{ y P 2 



(63) 



i H6 Transverse Y gradient Coil 

17 

1 8 Considering the current density symmetries amongst the planes of the self-shielded biplanar 

19 gradient coil design, as well as that the mutual energy is evaluated over two gradient coil 

20 representing the same axis, equation (63) has the form: 



19 



f ldadfi ^ 2 + fi2 cos{lk 0 )j;{a,p)j;\a 9 p) 
2 An -t P 

3 

4 The optimum phased array configuration occurs at the zo points where the value of the mutual 

5 inductance between the two sets is minimized. 
6 

7 X and Z gradient Coils 

8 Considering the current density symmetries amongst the planes of the self-shielded biplanar 

9 gradient coil design, as well as that the mutual energy is evaluated over two gradient coil 
fA representing the same axis, equation (63) has the form: 

it 9 

m 
» 

'■■'4 I 

B W M ^ =J^\ + ]dadp ^ a2 */ 2 cos{/3z 0 )J x a {a,P)j;\a,P) 
\i 4* J i P 



n 

ru 



iiy 

il[4 The optimum phased array configuration occurs at the z 0 points where the value of the mutual 

% inductance between the two sets is minimized. 

ny 

16 

17 Design and Results 

18 

19 l)Non-Shielded Biplanar Gradient Coils 

20 

21 i) Transverse Y biplanar gradient coil: 

22 For the non-shielded transverse biplanar gradient coils, where the gap between the two 

23 planes is chosen to be 2a=288mm, the constraint points for the minimization algorithm are 

24 shown in Table 1 . The gradient field at the center of the structure is demanded to be 40 mT/m 

25 with 10% on axis linearity at a distance of + 130mm from the center of the gradient and 20% of f 

26 axis uniformity inside a 25 cm Diameter Spherical Volume (DSV). 

20 



1 



2 
3 



4 
5 



a 



a 

m 

fi 
II 



Constraints for a Non-Shielded Biplanar Y gradient coil 


Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesla 


0.000 


0.001 


0.000 


0.0000400 


0.000 


0.130 


0.000 


0.0046800 


0.125 


0.001 


0.000 


0.0000320 


0.000 


0.001 


0.125 


0.0000320 


Table 1 



The application of ttfte Lagrange minimization technique to the set of constraints of the 
Table 1, generates a continuous current distribution for primary coil of the Transverse Non- 
Shielded Y biplanar coil. With tn^assist of the stream function technique, the continuous current 
densities for the two coils can be approximated by a 10 discrete loops with a common current of 
287.6 Amps. Table 2 illustrates all tire vital characteristics needed for the engineering and 
manufacturing phase of such a gradient coil. In addition figure 2 shows the on axis field behavior 
( on-axis linearity is represented by Figure 2a, off-axis uniformity is represented by Figure 2b )of 
the coil as it has been evaluated by applying the Biot-Savart law to the coil's discrete current 
pattern. 



Non-Shielded Y Biplanar Gradient Coil 


Property 




Gap 


288 mm 


Gradient Strength at mT/m/lOOA 


13.9 mT/m/lOOA 


Gradient Strength at 500 A 


69.5 mT/m 


450 A 


62.55 mT/m 


330 A 


45.87 mT/m 


Number of Discrete Loops 


10 


Cu thick, between adjacent loops 


> 12.5mm 


Total Inductance (cable included) 


95 uH 


Total Resistance(Cable included) 


55 mQ 



m 
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Rise Time at 2000V/500A 


24 usee 


700V/450A 


70 usee 




82 usee 


Slew Rate at 2000V/500A 


2865 T/m/sec 


700V/450 A 


880 T/m/sec 


400V/330A 


550 T/m/sec 


Duty Cycle (RMS Gradient) 




(500 A ref. point) at 4500 kW 


72% of the peak gradient(49.68mT/m) 


at 2000 kW 


49% of the peak gradient(33.81mT/m) 


at 1000 kW 


34% of the peak gradient(23.46mT/m) 



1 Table 2 



!; J ii) Transverse X biplanar gradient coil: 

□ 

For the non-shielded transverse biplanar gradient coils, where the gap between the two 
!;| planes is chosen to be 2a=284mm, the constraint points for the minimization algorithm are 
M| shown in Table 3. The gradient field at the center of the structure is demanded to be 40 mT/m 
B 7 with 10% on axis linearity at a distance of + 130mm from the center of the gradient and 20% off 
axis uniformity inside a 25 cm Diameter Spherical Volume (DSV). 



Constraints for a Non-Shielded Biplanar X gradient coil 


Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesla 


0.001 


0.000 


0.000 


0.0000400 


0.130 


0.000 


0.000 


0.0046800 


0.001 


0.125 


0.000 


0.0000320 


0.001 


0.000 


0.125 


0.0000320 



10 Table 3 

11 

The application of the Lagrange minimization technique to the set of constraints of the 

Table 3, generates a continuous current distribution for primary coil of the Transverse Non- 
Shielded X biplanar coil. With the assist of the stream function technique, the continuous current 
densities for the two coils can be approximated by a 10 discrete loops with a common current of 
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282.9 Amps). T?(ble 4 illustrates all the vital characteristics needed for the engineering and 
manufacturing phase of such a gradient coil. In addition figure 3 shows the on axis field behavior 
of the coil as it has ^een evaluated by applying the Biot-Savart law to the coil's discrete current 



pattern. 



Non-Shielded X Biplanar Gradient Coil 


Property 




Gap 


284 mm 


Gradient Strength at mT/m/lOOA 


15.62 mT/m/lOOA 


Gradient Strength at 500 A 
450 A 
330 A 


78.1 mT/m 

70.2 mT/m 
51.5 mT/m 


Number of Discrete Loops 


10 


Cu thick, between adjacent loops 


> 12.5mm 


Total Inductance (cable included) 


177uH 


Total Resistance(Cable included) 


59 mfi 


Rise Time at 2000V/500A 
700V/450A 


45 usee 
118 usee 
1 54 usee 


Linear Slew Rate at 2000V/500A 
700V/450 A 
400V/330A 


1738 T/m/sec 
593 T/m/sec 
335 T/m/sec 


Duty Cycle (RMS Gradient) 
(500 A ref. point) at 4500 kW 
at 2000 kW 
at 1000 kW 


68% of the peak gradient(53.1 lmT/m) 
45% of the peak gradient(35.37mT/m) 
32% of the peak gradient(25.0mT/m) 



Table 4 



Hi) Axial Z biplanar gradient coil: 

For the non-shielded Axial biplanar gradient coils, where the gap between the two planes 
is chosen to be 2a=292mm, the constraint points for the minimization algorithm are shown in 
Table 5. The gradient field at the center of the structure is demanded to be 40 mT/m with 10% on 
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1 

2 
3 



axis linearity at a distance of + 130mm from the center of the gradient and 20% of f axis 
uniformity inside a 25 cm Diameter Spherical Volume (DSV). 



Constraints for a Non-Shielded Biplanar Z gradient coil 


Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesla 


0.000 


0.000 


0.001 


0.0000400 


0.000 


0.000 


0.130 


0.0046800 


0.125 


0.000 


0.001 


0.0000325 


0.000 


0.125 


0.001 


0.0000325 



4 
5 



Table 5 



The application of th\ Lagrange minimization technique to the set of constraints of the 
Table 5, generates a continuous current distribution for primary coil of the Axial Non-Shielded Z 
biplanar coil. With the assist of the stream function technique, the continuous current densities 
for the two coils can be approximated by a 10 discrete loops with a common current of 298.15 
Amps. Table 6 illustrates all the vitakcharacteristics needed for the engineering and 
manufacturing phase of such a gradient coil. In addition figure 4 shows the on axis field behavior 
of the coil as it has been evaluated by applying the Biot-Savart law to the coil's discrete current 
li pattern. 
ll 
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Non-Shielded Z Biplanar Gradient Coil 


Property 




Gap 


292 mm 


Gradient Strength at mT/m/lOOA 


13.41 mT/m/lOOA 


Gradient Strength at 500 A 


67.08 mT/m 


450 A 


60.37 mT/m 


330 A 


44.27 mT/m 


Number of Discrete Loops 


10 


Cu thick, between adjacent loops 


> 12.5mm 


Total Inductance (cable included) 


73uH 
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Total Resistance(Cable included) 


61 mQ 


Rise Time at 2000V/500A 


19 usee 


700V/450A 


49 p.sec 


400V/330A 


64 usee 


Linear Slew Rate at 2000V/500A 


3530 T/m/sec 


700V/450 A 


1235 T/m/sec 


400V/330A 


698 T/m/sec 


Duty Cycle (RMS Gradient) 




(500 A ref. point) at 4500 kW 


61% of the peak gradient(40.66mT/m) 


at 2000 kW 


40% of the peak gradient(27.1mT/m) 


at 1000 kW 


29% of the peak gradient(19.2mT/m) 



Table 6 

2) Shielded Biplanar Gradient Coils 

i)Transverse Y biplanar gradient coil : 

For the shielded transverse biplanar gradient coils, the gap between the two primary 
coil's planes is chosen to be 2a=288mm, while the gap between the shielding planes is 
2b=337mm. The constraint points for the minimization algorithm are shown in Table 7. The 
gradient field at the center of the structure is demanded to be 40 mT/m with 10% on axis 
linearity at a distance of + 130mm from the center of the gradient and 20% of f axis uniformity 
inside a 25 cm Diameter Spherical Volume (DSV). 



Constraints for a Shielded Biplanar Y gradient coil 


Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesla 


0.000 


0.001 


0.000 


0.0000400 


0.000 


0.130 


0.000 


0.0046800 


0.125 


0.001 


0.000 


0.0000320 


0.000 


0.001 


0.125 


0.0000320 



Table 7 



The application of the Lagrange minimization technique to the set of constraints of the 
Table 7, generates a continuous current distribution for primary coil of the Transverse Shielded 



sist of th\st: 



14 Y biplanar coil. With the assist of th\ stream function technique, the continuous current densities 
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7 



1 for the primary cfoil can be approximated by a 15 discrete loops with a common current of 425.85 

2 Amps(figure 5). For a shielded Y biplanar gradient coil with parabolic returns, the current 

3 pattern for the primary coil is shown in Figure 6. while for the secondary coil its continuous 
current density can ba approximated by 10 loops (figure 7b)with the same constant current as the 
primary coils. For a smelded Y biplanar gradient coil with parabolic returns, the current 
pattern for the secondary coil is shown in Figure 7b. Table 8 illustrates all the vital 

7 characteristics needed for me engineering and manufacturing phase of such a gradient coil. In 

8 addition figure 8 shows the cm axis field behavior of the coil as it has been evaluated by applying 

9 the Biot-Savart law to the coil V discrete current pattern. - , , 
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Shielded Y Biplanar Gradient Coil 


Property 




Gap primary/shielding 


288 mm/3 3 7mm 


Gradient Strength at mT/m/lOOA 


9.39 mT/m/lOOA 


Gradient Strength at 500 A 


46.95 mT/ra 


450 A 


42.26 mT/m 


330 A 


3 1 .00 mT/m 


Number of Discrete Loops 


15/10 


Cu thick, between adjacent loops 


11. 4mm/ 15. 8mm 


Total Inductance (cable included) 


105uH 


Total Resistance(Cable included) 


95 mQ 


Rise Time at 2000V/500A 


27 usee 


700V/450A 


72 usee 


400V/330A 


94 usee 


Slew Rate at 2000V/500A 


1746 T/m/sec 


700V/450 A 


587 T/m/sec 


400V/330A 


330 T/m/sec 


Duty Cycle (RMS Gradient) 




(500 A ref. point) at 4500 kW 


47% of the peak gradient(22.27mT/m) 


at 2000 kW 


32% of the peak gradient(14.85mT/m) 
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at 1000 kW 


22% of the peak gradient(10.50mT/m) 


% Eddy Currents at 25cm DSV 


0.32% 



1 Table 8 

2 

3 iv) Transverse X biplanar gradient coil: 

4 For the shielded transverse biplanar gradient coils, where the gap between the two 

5 primary planes is chosen to be 2a=284mm, while the gap between the shielding planes is 

6 2b=341mm. The constraint points for the minimization algorithm are shown in Table 9. The 

7 gradient field at the center of the structure is demanded to be 40 mT/m with 10% on axis 

8 linearity at a distance of + 130mm from the center of the gradient and 20% off axis uniformity 

9 inside a 25 cm Diameter Spherical Volume (DSV). 

I J! 



Constraints for a Shielded Biplanar X gradient coil 


Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesla 


0.001 


0.000 


0.000 


0.0000400 


0.130 


0.000 


0.000 


0.0046800 


0.001 


0.125 


0.000 


0.0000320 


0.001 


0.000 


0.125 


0.0000320 



11 Table 9 

1L_ — ^ 

13 The application of the Lagrange minimization technique to the set of constraints of the 

14 Table , generates a continuous current distribution for primary coil of the Transverse Shielded X 

15 biplanar coil. With the assist of the stream function technique, the continuous current density for 

16 the primary coil can be approximated b\a 16 discrete loops with a common current of 339.2 

1 7 Amps(figure 9). The continuous current density of the shielding coil can be also approximated 

18 by a set of 1 1 discrete loops carrying the same current as the primary coil(figure 10). Table 10 

19 illustrates all the vital characteristics needed f^r the engineering and manufacturing phase of 

20 such a gradient coil. \ 
21 

22 
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Shielded X Biplanar Gradient Coil 


Property 




Gap primary/shielding 


284 mm/34 1mm 


Gradient Strength at mT/m/lOOA 


11.79mT/m/100A 


Gradient Strength at 500 A 


59.00 mT/m 


450 A 


53.10 mT/m 


330 A 


38.93 mT/m 


Number of Discrete Loops 


16/11 


Cu thick, between adjacent loops 


10.8mm/15.8mm 


Total Inductance (cable included) 


203 uH 


Total Resistance(Cable included) 


103 mQ 


Rise Time at 2000V/500A 


53 psec 


700V/450A 


140 psec 


400V/330A 


183 psec 


Slew Rate at 2000V/500A 


1132T/m/sec 


700V/450 A 


422 T/m/sec 


400V/330A 


213 T/m/sec 


Duty Cycle (RMS Gradient) 




(500 A ref. point) at 4500 kW 


45% of the peak gradient(26.55mT/m) 


at 2000 kW 


30% of the peak gradient(17.59mT/m) 


at 1000 kW 


21% of the peak gradient(12.44mT/m) 


% Eddy Currents at 25cm DSV 


0.37% 



Table 10 



v) Axial Z biplanar gradient coil: 

For the shielded Axial biplanar gradient coils, where the gap between the two primary 
planes is chosen to be 2a=292mm. The distance between the two shielding planes is 2b=345mm. 
The constraint points for the minimization algorithm are shown in Table 1 1 . The gradient field at 
the center of the structure is demanded to be 40 mT/m with 10% on axis linearity at a distance of 
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1 + 130mm from the center of the gradient and 7.5% of f axis uniformity inside a 25 cm Diameter 

2 Spherical Volume (DSV). 
3 



Constraints for a Shielded Biplanar Z gradient coil 


Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesla 


0.000 


0.000 


0.001 


0.0000400 


0.000 


0.000 


0.130 


0.0046800 


0.125 


0.000 


0.001 


0.0000394 


0.000 


0.125 


0.001 


0.0000394 



4 Table 11 

5 



; (5 The application of the Mgrange minimization technique to the set of constraints of the 

Table 11, generates a continuou^current distribution for primary coil of the Axial Shielded Z 
biplanar coil. With the assist of th^stream function technique, the continuous current densities 
! j for the primary coil can be approximated by 16 discrete loops with a common current of 420.9 
jib Amps(figure 1 1). the continuous curreVt density for the shielding coil can be approximated by 12 
;il 1 discrete loops (figure 12) Table 12 illustrates all the vital characteristics needed for the 
jijf2 engineering and manufacturing phase of s^ch a gradien t coil. 

nil 



Shielded Z Biplanar Gradient Coil 


Property 




Gap primary/shielding 


292 mm/345mm 


Gradient Strength at mT/m/lOOA 


9.50mT/m/100A 


Gradient Strength at 500 A 


47.51 mT/m 


450 A 


42.76 mT/m 


330 A 


31.36 mT/m 


Number of Discrete Positive Loops 


16/12 


Cu thick, between adjacent loops 


11.0mm/l 5.8mm 


Total Inductance (cable included) 


108 uH 
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Total Resistance(Cable included) 


107 mQ 


Rise Time at 2000V/500A 
700V/450A 
400V/330A 


28 usee 
75 usee 
98 usee 


Slew Rate at 2000V/500A 
700V/450 A 
400V/330A 


1712T/m/sec 
573 T/m/sec 
320 T/m/sec 


Duty Cycle (RMS Gradient) 
(500 A ref. point) at 4500 kW 
at 2000 kW 
at 1000 kW 


44% of the peak gradient(20.68mT/m) 
29% of the peak gradient(13.79mT/m) 
21% of the peak gradient(9.75mT/m) 


% Eddy Currents at 25cm DSV 


0.33% 



g Table 12 



3" 3) Non-Shielded Uniplanar Gradient Coils 

m 

; : 

5 i) Transverse X, Y uniplanar gradient coil 

For the non-shielded transverse uniplanar gradient coils, where the gap between the plane 

i y 

ITU from the geometric center of the structure of the gradient coil is chosen to be -a'=-6mm, the 

ru 

;8j constraint points for the minimization algorithm are shown in Table 13. The gradient field at the 

^ center of the structure is demanded to be 40 mT/m with 10% on axis linearity at a distance of + 

10 100mm from the center of the gradient and 20% of f axis uniformity inside a 25 cm Diameter 

1 1 Spherical Volume (DSV). 
12 



Constraints for a Non-Shielded Uniplanar transverse(X,Y) gradient coils 


Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesla 


0.000 


0.001 


0.100 


0.0000400 


0.000 


0.060 


0.100 


0.0024000 


0.000 


0.200 


0.100 


0.0072800 


0.180 


0.001 


0.100 


0.0000360 
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0.000 


0.001 


0.150 


0.0000410 


0.000 


0.001 


0.200 


0.0000350 



Table 13 



The application of the Lagrange minimization technique to the set of constraints of the 
Table 13, generates a continuous current distribution for primary coil of the Transverse Non- 
Shielded Y uniplanar coil. With the assist of the stream function technique, the continuous 
current densities for the coil can be approximated by a 10 discrete loops with a common current 
of 41 1.2 Amps (figure 13). Table 14 illustrates all the vital characteristics needed for the 
engineering and manufacturing phase of sucl\a gradient coil. In addition figure 14 shows the on 
axis field behavior of the coil as it has been evaluated by applying the Biot-Savart law to the 
coil's discrete current pattern. \ 



Non-Shielded Transverse (X,Y) Uniplanar Gradient Coil 


Property 








Gradient Strength at mT/m/lOOA 


9.73 mT/m/lOOA 


Gradient Strength at 500 A 
450 A 
330 A 


48.64 mT/m 
43.77 mT/m 
32.11 mT/m 


Number of Discrete Loops 


10 


Cu thick, between adjacent loops 


> 10.2mm 


Total Inductance (cable included) 


295^H 


Total Resistance(Cable included) 


73 mfi 


Rise Time at 2000V/500A 
700V/450A 
400V/330A 


75 usee 
199 usee 
259 (asec 


Slew Rate at 2000V/500A 
700V/450 A 
400V/330A 


648 T/m/sec 
220 T/m/sec 
124 T/m/sec 
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Duty Cycle (RMS Gradient) 




(500 A ref. point) at 4500 kW 


53.4% of the peak gradient(26 mT/m) 


at 2000 kW 


36% of the peak gradient(17.33mT/m) 


at 1000 kW 


25% of the peak gradient(12.26mT/m) 


Table 14 



1 

2 
3 
4 
5 
6 
7 
8 

$>] inside a 20 cm Diameter Spherical Volume (DSV). 



ii) Axial Z Uniplanar gradient coil: 

For the non-shielded Axial Uniplanar gradient coil, the gap between the geometric center 
of the structure of the gradient coil and the planar surface is chosen to be -a--6mm. The 
constraint points for the minimization algorithm are shown in Table 15. The gradient field at the 
center of the structure is demanded to be 40 mT/m with 20% on axis linearity at a distance of + 
200mm from the center of the gradient with no-rollover point and 20% of f axis uniformity 
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14 



16 
17 
18 
19 



Constraints for a Non-Shielded Uniplanar Z gradient coil 


Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesla 


0.000 


0.000 


0.031 


0.0000400 


0.000 


0.000 


0.231 


0.0062800 


0.180 


0.000 


0.031 


0.00003600 


0.000 


0.200 


0.031 


0.00003200 



Table 15 



The application of the Lagrange minimization technique to the set of constraints of the 
Table 15 generates a continuous current distribution for primary coil of the Axial Non-Shielded 
Z uniplanar coil. With the assist of the stream function technique, the continuous current 
densities for the two coils can be approximated by a 10 discrete loops with a common current of 
316.1 Amps(figure 15). Table 16 illustrates all the vital characteristics needed for the 
engineering and manufacturing phase of such a\gadient coil. ^ 
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Non-Shielded Z Uniplanar Gradient Coil 


Property 








Gradient Strength at mT/m/lOOA 


12.65 mT/m/lOOA 


Gradient Strength at 500 A 
450 A 
330 A 


63.25 mT/m 
56.93 mT/m 
41.75 mT/m 


Number of Discrete Loops 


13 


Cu thick, between adjacent loops 


> 12.5mm 


Total Inductance (cable included) 


241uH 


Total Resistance(Cable included) 


57 mQ 


Rise Time at 2000V/500A 
700V/450A 
400V/330A 


62 usee 
160 psec 
209 psec 


Linear Slew Rate at 2000V/500A 
700V/450 A 
400V/330A 


1034T/m/sec 
354 T/m/sec 
200 T/m/sec 


Duty Cycle (RMS Gradient) 
(500 A ref. point) at 4500 kW 
at 2000 kW 
at 1000 kW 


63% of the peak gradient(40.0mT/m) 
42% of the peak gradient(26.7mT/m) 
30% of the peak gradient(18.9mT/m) 



1 Table 16 

2 4) Shielded Uniplanar Gradient Coils 

3 i)Transverse (X,Y) uniplanar gradient coil with parabolic cut-offs: 

4 For the shielded transverse biplanar gradient coils, the gap between the geometric center 

5 of the entire gradient structure and the planar surface of the primary coil is chosen to be -a'=- 

6 6mm. The gap between the geometric center of the structure and the planar surface of the 

7 shielding coil is chosen to be -b'=-70mm. The constraint points for the minimization algorithm 

8 are shown in Table 1 7. The gradient field at the center of the structure is demanded to be 40 

9 mT/m with 20% on axis linearity at a distance of 
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+ 200mm from the center of the gradient. In addition, the off-axis uniformity of the gradient field 
is restricted to within 20% from the field's ideal value inside a 25 cm Diameter Spherical 
Volume (DSV). 



Constraints for a Shielded Uniplanar transverse(X,Y) gradient coils 


Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesla 


0.000 


0.001 


0.075 


0.0000400 


0.000 


0.100 


0.075 


0.0040000 


0.000 


0.300 


0.075 


0.0140400 


0.180 


0.001 


0.075 


0.0000380 


0.000 


0.001 


0.13 


0.0000390 


0.000 


0.001 


0.180 


0.0000320 



Table 17 



The application of the Lagraritee minimization technique to the set of constraints of the 
Table 17 generates a continuous current distribution for primary coil of the Transverse Shielded 
Y biplanar coil. With the assist of the stream function technique, the continuous current densities 
for the primary coil can be approximated b\ a 16 discrete loops with a common current of 324.17 
Amps(figure 16a), while for the secondary coil its continuous current density can be 
approximated by 8 loops (figure 16b)with the same constant current as the primary coils. Table 
18 illustrates all the vital characteristics needed iW the engineering and manufacturing phase of 
such a gradient coil. In addition figure 16c shows tf^e on axis field behavior of the coil as it has . 
been evaluated by applying the Biot-Savart law to thXcoil's discrete current pattern. 



Shielded Transverse (X,Y) Uniplanar Gradient Coil 


Property 








Gradient Strength at mT/m/lOOA 


12.33 mT/m/lOOA 


Gradient Strength at 500 A 
450 A 
330 A 


61.69 mT/m 
55.52 mT/m 

40.70 mT/m 
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Number of Discrete Loops 


16/8 


Cu thick, between adjacent loops 


12.4mm/ 15. 8mm 


Total Inductance (cable included) 


302^iH 


Total Resistance(Cable included) 


122 mQ 


Rise Time at 2000V/500A 
700V/450A 
400V/330A 


78 usee 
210 usee 
277 usee 


Slew Rate at 2000V/500A 
700V/450 A 
400V/330A 


792 T/m/sec 
264 T/m/sec 
147 T/m/sec 


Duty Cycle (RMS Gradient) 
(500 A ref. point) at 4500 kW 
at 2000 kW 
at 1000 kW 


40% of the peak gradient(24.95mT/m) 
27% of the peak gradient(16.64mT/m) 
19% of the peak gradient(l 1.76mT/m) 


% Eddy Currents at 25cm DSV 


0.39% 



* Table 18 

a 
a 



it 

[4j Hi) Axial Z biplanar gradient coil: 

V'a I 
: [ i~ 

For the shielded Axial Uniplanar gradient coils, the gap between the geometric center of 

6 the entire gradient structure and the planar surface of the primary coil is chosen to be -a'=-6mm. 

7 The gap between the geometric center of the structure and the planar surface of the shielding coil 

8 is chosen to be -b --70mm. The constraint points for the minimization algorithm are shown in 

9 Table 19. The gradient field at the center of the structure is demanded to be 40 mT/m with 20% 

10 on axis linearity at a distance of ± 130mm from the center of the gradient and 7.5% of f axis 

1 1 uniformity inside a 25 cm Diameter Spherical Volume (DSV). 
12 



Constraints for a Shielded Uniplanar Z gradient coil 


Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesla 


0.000 


0.000 


0.031 


0.0000400 
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0.000 


0.000 


0.231 


0.0062800 




0.180 


0.000 


0.031 


0.00003600 




0.000 


0.200 


0.031 


0.00003200 


1 


Table 19 








2 











3 The application of tj^e Lagrange minimization technique to the set of constraints of the 

4 Table 19, generates a continuous current distribution for primary coil of the Axial Shielded Z 

5 biplanar coil. With the assist of me stream function technique, the continuous current densities 

6 for the primary coil can be approximated by 21 discrete loops with a common current of 385.71 

7 Amps(figure 17a). The continuous current density for the shielding coil can be approximated by 

8 13 discrete loops (figure 17b) Table 20 illustrates all the vital characteristics needed for the 
; 9 engineering and manufacturing phase o\such a gradient coil. 
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Shielded Z Uniplanar Gradient Coil 


Property 








Gradient Strength at mT/m/lOOA 


10.37 mT/m/lOOA 


Gradient Strength at 500 A 
450 A 
330 A 


51.85 mT/m 
46.67 mT/m 
34.22 mT/m 


Number of Discrete Positive Loops 


21/13 


Cu thick, between adjacent loops 


9.4mm/ 15.8mm 


Total Inductance (cable included) 


306 uH 


Total Resistance(Cable included) 


125 mQ 


Rise Time at 2000V/500A 
700V/450A 
400V/330A 


80 usee 
213 usee 
281 usee 


Slew Rate at 2000V/500A 
700V/450 A 


656 T/m/sec 
218 T/m/sec 
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400V/330A 


122 T/m/sec 


Duty Cycle (RMS Gradient) 
(500 A ref. point) at 4500 kW 
at 2000 kW 
at 1000 kW 


40% of the peak gradient(20.97mT/m) 
27% of the peak gradient(13.98mT/m) 
19% of the peak gradient(9.89mT/m) 


% Eddy Currents at 25cm DSV 


0.19% 



1 Table 12 

2 

3 Mutual Inductance Calculations 

4 

5 Biplanar Gradients 

6 

lis!: 

% s The mutual energy between two biplanar gradients of the same axis was evaluated based 

£p on the equations (64) and (65). In particular for this example the mutual energy between two 

: :i : 
•: "5? 

iSg biplanar transverse Y gradients was evaluated. For the first biplanar set, the gap between its 

10; primary planes was set at 2a=288mm, while the gap between its shielding planes was set at 

IV 2b=341mm. For the second biplanar set, the gap between the primary planes was set at 

128 2a=282mm, while the gap between the two secondary planes is set at 2 b =292mm. The mutual 

1| energy of these two sets as a function of shifting one set with respect to the other along the z 

iM direction. The optimum (zero) mutual energy of these two sets occurs at the point zo=216mm 

□ 

1| from the geometric center of the fir gradient coil set. 
16 

17 Uniplanar Gradients 
18 

19 The mutual energy between two uniplanar gradients of the same axis was evaluated based 

20 on the equations (21). In particular for this example the mutual energy between two uniplanar 

21 axial Z gradients was evaluated. For the first uniplanar set, the gap on the primary plane was set 

22 at -a—6mm, while the gap of the shielding plane was set at -b --70mm. For the second biplanar 

23 set, the gap for the primary plane was set at -a--3mm, while the gap between the two secondary 

24 planes is set at - b --75mm. The optimum (zero) mutual inductance for these two set occurs at 

25 the point zo=245mm from the geometric center of the first gradient coil set. 
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2 The present invention describes a method of designing and manufacturing an Open 

3 Architecture gradient coil set which is suitable for magnetic resonance systems with vertically 

4 directed or horizontally directed fields. The concept of designing of an open architecture gradient 

5 coil system offer a significant advantage over existing geometries since, the present gradient 

6 design provides superior performance in terms of gradient strength and field linearity and 

7 uniformity while the structure of the coil is an open architecture focusing of reducing patient 

8 claustrophobia. The invention can be implement in gradient coil designs for open magnet 

9 geometries with either vertically and or horizontally directed fields. In particular, for magnets 

10 with vertically directed fields, besides the open architecture design of the proposed gradient 

1 1 system a significant gradient increase in the strength of the local gradient field up to 55mT/m 

12 with slew rates exceeding 1000 T/m/sec in an imaging volume covering the size of the human 
IX head (25 cm DSV). Another significant advantage of the design is the ability to control and 

tM , minimize the dB/dt and eddy current levels for shielded design. In addition, the open architecture 

1 1 design can be paired in a phased array configuration with alternative shielded or not open 

architecture gradient coil set in such a configuration that covers the human body head to toe. The 

13 significant advantage of such configuration is the significant reduction of dB/dt effects and the 
1§ utilization of high strength non-claustrophobic gradient coils for high speed imaging of the 

l| human body with excellent gradient field uniformity and linearity characteristics in the local 

20 predetermined gradient volume. For open magnets with horizontally directed fields (double- 

: : sat 

2!j donut designs), the invention provides a competitive advantage over existing gradient designs 

22 since it can deliver gradient fields with strengths two to three fold grater than already existing 

23 designs and the ability to cover the entire human body in a phased array configuration with 

24 minimized dB/dt effects and induced eddy currents on the magnet's shields or metallic surfaces 

25 during the localization of the gradient field. The open architecture can be modified with 

26 parabolic cut-offs in order to incorporate the human shoulders for the portion of the open 

27 architecture gradient system which is suitable for imaging the human head. Furthermore, the 

28 design can be incorporated as an insertable configuration and can be fully integrated with any 

29 volumetric or surface phased array or singular coil configuration for imaging portions or the 

30 entire human body. In another aspect of the invention the design can be secured into the structure 

31 of its corresponding design. Although the design incorporates a limited curvature for the 
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1 perspective gradient coil configurations, the curvature can be easily modified for those skilled to 

2 the art to any unspecified shape suitable for the any targeted application. In addition, for a self- 

3 shielded design the side planes of the secondary (shielded) coil can be mechanically shifted and 

4 can be adjusted with respect the primary planes of the coil in an array configuration with zero 

5 mutual inductance in such a way that the imaging volume of the coil is extended. 

6 The present invention provides a gradient coil set that minimizes patient anxiety while 

7 delivering high peak, uniformity and slew rate gradient coils for Large imaging volumes. The 

8 combination of both biplanar and uniplanar gradients achieve a gradient magnetic field with 

9 superior field qualities. The coil set may have the side planes with parabolic cut offs in order to 

10 assist in patient comfort for applications that are focused on imaging the head, the neck, the 

1 1 spine, the heart and the upper torso area. Due to the phased array design low dB/dt levels can be 

12 achieved. In the gradient coil set and especially for the shielded designs, the outer shielding 

ii planes can be moved to form a phased array configuration with their perspective primary planes 

l| corresponding to the same gradient coil axis. The movement can be done via mechanical means 

1:5 or any other means. An unlimited spectrum of FMRI and Interventional applications can be 

""i 

16 performed with the he "Open Face Architecture" gradient coil, because of its openness and the 

ll high peak and slew rate values for the gradient field. 

W Referring to FIG. 1, a gradient coil set 10 includes a uniplanar Z-gradient coil 12 and a 

biplanar X-gradient coil 14 and a biplanar Y-gradient coil 1 6. Each of coils may include a 

2® respective shield coil 18, 20, 22. The coil set 10 has an open Z-axis face opposite the coil 12. 

Q 

M The biplanar coils may be provided with shoulder reliefs 24. The structure of the biplanar coils 

22 14, 16 may be conjoined with the structure of the uniplanar coil 12, for example, as shown. 




Referring to FIG. 18, ^gradient coil set 30 includes self-shielded gradient coils 32, 34, 36 
and a RF head coil 38 

Referring to FIG. 19, a gradient coil set 40 includes un-shielded gradient coils 42, 44, 46 
and a RF head coil 48. 

Referring to FIG. 20, a gradient coil set 50 includes self-shielded gradient coils 52, 54, 
56, a RF head coil 58 and a CTL spine ®^il 59 integrated into the structure of the gradient coil set 
50. 

Referring to FIG. 21, a gradient coil ^et 60 includes un-shielded gradient coils 62, 64, 66, 
a RF head coil 68 and a CTL spine coil 69 integrated into the structure of the gradient coil set 60. 
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1 Referring tfi) FIG. 22, a gradient coil set 70 includes Z-gradient uniplanar coils 72, 74 and 

2 X-Y biplanar gradient coils 76, 78 in a phased array configuration. A head RF coil 79 may be 

3 included. 
y \ 

Referring to FI(1 23, an open Z-axis face coil set 90 includes a knee and upper thigh RF 

ty» 5 coil 92 integrated therewn* 

6 Referring to FIG. 24^ an open Z-axis face coil set 100 includes a knee and foot RF coil 

7 102 integrated therewith. 



8 It should be evident that this disclosure is by way of example and that various changes 

9 may be made by adding, modifying or eliminating details without departing from the fair scope 

10 of the teaching contained in this disclosure. The invention is therefore not limited to particular 

1 1 details of this disclosure except to the extent that the following claims are necessarily so limited. 
12 

If 4 
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